Magneto-spin Hall conductivity of a two-dimensional electron gas 
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It is shown that the interplay of long-range disorder and in-plane magnetic field gives rise to an 
out-of-plane spin polarization and a finite spin Hall conductivity of the two-dimensional electron gas 
in the presence of Rashba spin-orbit coupling. A key aspect is provided by the electric-field induced 
in-plane spin polarization. Our results are obtained first in the clean limit where the spin-orbit 
splitting is much larger than the disorder broadening of the energy levels via the diagrammatic 
evaluation of the Kubo-formula. Then the results are shown to hold in the full range of the disorder 
parameter appr by means of the quasiclassical Green function technique. 



It is well established that the peculiar linear-in- 
momentum dependence of the Rashba (and of Dres- 
selhaus) spin-orbit coupling leads in a two-dimensional 
electron gas (2DEG) to the vanishing of the spin Hall 
conductivity [TJ [21 [31 2]. This can be directly recognized 
by considering the continuity-like equation for the in- 
plane spin polarization, where the spin-nonconserving 
terms can be written as the spin current associated to 
the out-of-plane spin polarization and to the spin Hall 
effect [31 [3- In this paper, we show that the inter- 
play of an in-plane magnetic field, B, taken parallel to 
the electric field, £, (say along the e^, axis) and long- 
range disorder changes this behavior providing then a 
potential handle on the spin Hall effect. In particular, 
we show that while the out-of-plane spin polarization is 
linear in the magnetic field, the spin Hall conductivity 
is quadratic. Our analysis is valid in the standard good 
metallic regime cft/K^ 1 with spin-orbit effects taken 
into account to first order in a/vp- Here a and r are 
the spin-orbit coupling and the elastic quasiparticle life- 
time due to impurity scattering, respectively, while vp, 
Pf and ep — vppp/2 are the parameters of the 2DEG in 
the absence of spin-orbit coupling. 

Our proposal of a magnetic field-induced spin Hall ef- 
fect differs from related previous suggestions both for 
the analytical treatment of it [8 and for the microscopic 
mechanism responsible of the effect [5]. The difference 
of our proposal with respect to Ref. [5] is closely related 
to the electric field-induced in-plane spin polarization, 
which for short-range disorder scattering is given by [10' 

Sy = ~NQa\e\£xT. (1) 

In Eq.([T]) A^o — »Ti/(27r) is the free density of states of the 
2DEG in the absence of spin-orbit interaction. Contrary 
to what one could expect, the generalization of Eq.([lJ for 
long-range disorder, as we will show later, is not the re- 
placement of the elastic quasiparticle lifetime r with the 
transport time, Ttr, as it has been assumed in Ref. [9]. We 
find then that long-range disorder leads to a non-trivial 
modification of the effective Bloch equations and even- 



tually yields out-of-plane spin polarization and spin Hall 
effect. In contrast to Ref. |H] we do not have to assume 
a non-parabolicity of the energy bands or an energy de- 
pendence of the scattering probability. 

Our analysis is carried out in two steps. In the first 
step, we calculate the out-of-plane spin polarization using 
the diagrammatic approach of Ref. [3J , valid in the clean 
limit when the spin-orbit splitting is much larger than 
the disorder-induced broadening, 2app ^ t~^. In order 
to make contact with the analysis of Ref. [H] performed in 
the opposite dirty limit {app <^ t^^), we present, in the 
second step, a derivation based on the Eilenberger equa- 
tion for the quasiclassical Green function in the presence 
of spin-orbit coupling [11]. The advantage of so doing is 
that the analysis is valid for an arbitrary value of the pa- 
rameter 2appT and also allows to determine the effective 
Bloch equations for the spin density. 

The Hamiltonian of a 2DEG perpendicular to the e^- 
axis reads 

ff= |^+b(p)-<T + T/(x) (2) 

where b(p) = ap x — lUsGx is the effective magnetic 
field including both the Rashba spin-orbit coupling and 
the external magnetic field. In Eq.([2|, V{x) describes the 
potential scattering from the impurities and, according to 
the established procedure in the literature, will be taken 
as a random variable. 

The possibility of a non- vanishing spin Hall conductiv- 
ity in the presence of an in-plane field may be appreciated 
by considering the equation of motion for the ej,-axis (in- 
plane) spin polarization, which yields 

ds d 

^ + ^ • - -2majXy + 2uss,, (3) 

where jl ^ is the e-y-axis component of the spin current 
that is polarized along the e^-axis. Under stationary and 
uniform conditions, the above equation implies, in the 
absence of the magnetic field, a vanishing spin current 
and hence a vanishing spin Hall conductivity, with the 
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FIG. 1: Diagrams to be evaluated in zeroth and first order 
in the magnetic field. The dashed line ending with a cross 
indicates the magnetic field insertion. Under impurity average 
the (spin density and charge current) vertices are dressed by 
the standard ladder resummation. 



latter defined by j^y — (^s,h^x- In the presence of an 
in-plane magnetic field one may have a bulk spin Hall 
current determined by 



Js 



UJs 

am 



(4) 



In the following we evaluate to first order in the mag- 
netic field which implies that of y to second order. 

As anticipated, we begin by sketching the calculation 
performed with the diagrammatic approach. To linear 
order in the electric field, the Kubo formula for the zero- 
temperature expression of the out-of-plane spin polariza- 
tion reads 



^ J2 Tr. kG«(p)j-G^(p)l |e|£,, (5) 



27r 



where the Green functions can be obtained from Eq.(|2|, 
the vertices are Sz — {l/2)az, = (Pk/'7i)o'o ~ aay, and 
the trace is over the associated spin indices. In Eq.(|5]), 
the bar indicates the average over the impurity potential, 
which, at the level of the self-consistent Born approxima- 
tion, yields the self-energy 



p' 



(6) 



|F(q)p being the Fourier transform of y(x)V(x'). In 
order to consider the effect of long-range disorder, we 
expand the above scattering probability as 



\V\^ = Vo + 2ViCos{ip-ip') + 2V2Cos{2ip-2ip') + - 



(7) 



where ip—(p' is the angle between the two momenta p and 
p'. The harmonics, Vq, Vi, V2 are functions of |p| and 
|p'|. In the following we will ignore this dependence [12] 
and take both momenta at pp. To first order in the 
magnetic field the diagrams to be evaluated are shown in 
FigjT] Notice that all the vertices and propagators ap- 
pearing in the diagrams must be evaluated at zero mag- 
netic field. It is then convenient, following Ref.^, to use 
the disorder-free Hamiltonian eigenstates 



|P±> 



1 



(±iexp(-i(p)|pT) + |pi)) 



(8) 



corresponding to the eigenvalues E± = p^/2m± ap with 
taii(ip) = py/px- In terms of the transformation matrix, 
J7, defined by Eq.(|8|, the Pauli matrices transform as 



UayW = 
Ua,W = 



Py 



PxO-z +Py<7y 
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cos{(p) 
sin((p) 



(9) 



As a consequence the spin density vertex, the magnetic 
field insertion and the charge current vertex become 
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(10) 

(11) 
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Upon impurity averaging the spin and charge vertices get 
renormalized. In terms of the renormalized quantities, 
Eq.(|5]) becomes to first order in the Zeeman field 



\e\£x 
in 



^ ] Px^J■i^ fj.fj.G Tpfj^Gp )G ^ Jc,fifn 
p 

(13) 

where the first (second) term in the brakets refers to the 
magnetic field insertion in the top (bottom) Green func- 
tion line and fJ. — i: labels the eigenstates. The quantities 
r^^/ and Jcfifj.' are the dressed vertices corresponding to 
Sz and j^, respectively, and the Green functions are eval- 
uated via the self-energy given in Eq.([6]). Apart from 
the spin vertex T^^/, all the other quantities have been 
evaluated in Ref. [3 , where it has been shown that the 
off-diagonal matrix element Jc.^/z vanishes and the self- 
energy is diagonal in the eigenstate basis with 

(14) 

The spin vertex obeys the equation 

p' 



which yields 



iM Vi 
2apFT Vq 



(15) 



(16) 



By using G^G^ = iT^(G^ — G^), integrating over the 
energy, ^ — p^ /2m — ji, and keeping terms up to order 
a/vp, Eq.(13l becomes 



\e\£x 



UJs 



8ah 



where N± = Nq{1 q= a/vp), P± = Pf{^ T oi/vp) are 
the density of states and the Fermi momentum of the 
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two spin subbands and J±±{p±) = J±'Px are the Fermi- 
surface expressions of the charge vertices. Finally, bor- 
rowing from Ref. [3] the expression for the vertices 



J± = Vp 



Vo 



a V0 + V2 



VF Vo - V2 

one gets the out-of-planc spin polarization 



1 UJs N0 V1-V2 

2 apF pp Vq - V2 



and, via Eq.(|4|, the spin Hall conductivity 

V2 



^3. 

in 



apF , 



V0-V2 



(18) 



(19) 



(20) 



Remarkably the above central result shows that the sign 
of the spin Hall conductivity depends on the relative 
strength of the harmonics of the scattering probability. 
This may explain the sign change in the numerical evalua- 



tion of Ref . [8| . On the other hand, Eq.(19l is inconsistent 
with Ref. [9^ , where for the disorder model of Eq. ([t]) one 
would expect a vanishing out-of-planc polarization. 

Recall that we derived Eqs.(19l and pOj ) in the clean 
limit. We now rederive the above results by means of 
the kinetic equations approach of Ref. [TT] and find that 
Eqs.(19l and (20 1 are valid for all values of the disorder 
parameter appT. Furthermore we will derive the effective 
Bloch equations in the dirty limit. We start with the 
Eilenberger equation 



(21) 



for the quasiclassical Green function {g = gtjt2(p;x)) 

9=^1 deGt,t,(p,x), G = ^ (22) 

where Gtjf2(p,x) is the Wigner representation of the 
Green function, which has both matrix structure in the 
Keldysh (denoted by the check symbol) and spin spaces. 
[,] and {,} indicate commutator and anticommutator. 
As for the diagrammatic approach, the index /i = ± la- 
bels the two spin subbands (b± = b(p±)). In integra- 
tions like in Eq. ( 22 1 the corresponding poles in the Green 



functions yield the two-component decomposition of the 
quasiclassical Green function 



9± 



l{i(l±bo.^),^}. 



(23) 



The " 0" subscript denotes evaluation at the Fermi surface 
in the absence of spin-orbit coupling. In the following we 



are going to use Eq.(21l to first order in the parame- 
ter |bo|/eF- The connection to the physical observables 



is made by integrating over the energy e, which is the 
Fourier conjugated variable of the time difference ti — ^2- 
For instance, the out-of-plane spin density is given by the 
angular average of the Keldysh component [13] 



de(Tr(cr^5)), 



2tt 



(24) 



In order to solve, to linear order in the electric field, the 



Keldysh component of the Eilenberger equation (21), we 
use the minimal substitution 9x3 ~^ —\e\£x^xdegeq where 
= tanh(e/2r)(.gf,-54) with^ = = l-d^WcT 
is the equilibrium quasiclassical Green function. As in 
the diagrammatic treatment previously developed, we 
find it convenient to transform the equations to the eigen- 
state basis via Eq.([9|. After expressing the quasiclassical 
Green function as a four-dimensional column vector 



UgU^ 



30 53 51 92 



(25) 



the Eilenberger equation (211 can be then written as a 



linear system of four equations for the components of g 
dt~g = --{Mo + Mi)~g+-{l + N) {K~g) + So + Su (26) 

T T 



where, 
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(27) 



(28) 



(29) 



In Eq.(26l, (...) denotes angle integration over with 
the scattering kernel that can be expandend into angular 
harmonics as 

K{ip-if') = if(")+cos((^-^')^^°^+sin(</3-(/3')^^''^ + - • • 

(30) 

each coefficient being itself a matrix 



and 
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Finally the source electric-field dependent terms are 



So = E 



( Px \ 





, Si 



vfPf 



-E 



( \ 

PxPy 


V pI J 



with E = \e\£xVFd^{2ta.nh{e/2T)). Notice that, consis- 
tently with the accuracy we are working, one may use 
for the charge density component the solution obtained 
in the absence of both spin-orbit coupling and magnetic 
field 



jc,x ~ (Pxgo ) 



;;(0)\ _ 



1 Vn 



2Vo~Vi 



tE, 



where the characteristic transport time renormalization 
Ttr = tVo/{Vo — Vi) appears. We seek now a stationary 
solution of Eq. pS] ) which is evaluated in first order in the 
magnetic field g — g^^'' + g^^^ -|- . . . . We then get 



(.9(1)) = (Afo - if(°))-i(r(5i) - (Afi5(°))). 



(32) 



According to the transformation of Eq.(10l, the out-of- 
plane spin polarization is related to 
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vfPf 
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1 Vo-Vi 



apFT 



Vn 



{Pxgs 



(33) 

which is expressed in terms of (pxgs) evaluated at zero 
magnetic field. This latter quantity is nothing but the in- 
plane spin polarization (cf. the second line of Eq.([9|). By 



multiplying the second component of the system Eq. ( 26 ) 



by Px and performing the angle average, one obtains the 
generalization, for long-range disorder, of the Edelstein 
result ^lOj 



(Pxgl 



Vn 



vf Vq - V2 



-tE. 



(34) 



Finally, by using Eqs.( |33|34[ ) into Eq.(|24l one recovers 
the result ( 19 1 of the diagrammatic approach, which is 
now manifestly valid for any strength of the disorder. To 



understand the meaning of Eq. ( 34 1 , it is useful to recall 



the origin of the in-plane polarization jTO]: In the pres- 
ence of an electric field the Fermi surface is shifted by 
5px ~ \r\£xT- As a result the total spin of the electrons 
neither in the plus nor in the minus band adds up to zero. 
Although the contributions of both bands tend to cancel, 
a finite spin polarization remains due to the a/vF cor- 
rections in the density of states. For long-range disorder, 
the Fermi surface shift is proportional to the transport 
time Ttr, so one might expect the transport time also 
in the in-plane spin polarization. However due to the 
a/vp corrections each band has its own effective trans- 
port time, Ttr,± = J±t±/vf and the explicit result reads 
(so = aNf)\e\£xT) 



VF 



{N+Ttr,+ - N_Ttr,-)\e\£x 



Vo 



V0-V2 



so- (35) 



At last we study the combined effect of magnetic and 
electric field in the diffusive regime, uJsT,apFT ^ 1. 
The effective Bloch equations for the spin density are 

dtsx = -T-\sx-sl'^) (36) 
dtsy = [sy + soVo/{Vo - V2)-'] + 2u;ss, (37) 

dtSz = -2t^^Sz - 2u!s [sy + soTtr/r] , (38) 



where ^ — 2{apF)^Ttr and s^'^ = Ngajs. Eqs.(36l and 



( 38 1 agree with what was found in Ref . [9] , the only dif- 
ference is in the term proportional to the electric field 



(i.e. So) in Eq.(37). The stationary spin polarization as 



(31) a function of magnetic field is now determined as 



VqSq 14 
V0-V2 
Voso Vi 



KVo 
1 + 



2^2 



Vo 



(40) 



showing an out-out-plane contribution as observed ex- 
perimentally in Ref. |14l . 

In conclusion, we have shown that the combined effect 
of an in-plane magnetic field, long-range disorder and 
spin-orbit coupling gives rise to an out-of-plane spin po- 
larization and finite spin Hall conductivity, whose value 
does not depend on the concentration of defects as long 
as the 2DEG is in the metallic regime. To obtain the 
correct value of the electric-field induced in-plane spin 
polarization it is essential to take into account the differ- 
ent transport times in the two spin-orbit splitted bands. 
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